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Abstract 


The molten carbonate fuel cell (MCFC) is a promising device for stationary power and heat supply. In this work, the temperature control in 
MCFC stack is studied. The operating temperature of the stack is an important controlled variable, which impacts the performance of the MCFC 
due to thermal cycling. In order to improve the generating performance of the MCFC, prolong its life and guarantee safety, credibility and low 
degradation of the fuel cell, the temperature of MCFC must be controlled efficiently. Based on analyzing the temperature dynamic model of MCFC, 
a model-based variable structure control (VSC) controller is designed for MCFC in this paper. The optimal parameters of VSC are selected by 
genetic algorithm (GA). Numerical results demonstrate the effectiveness and advantages of this approach. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


The molten carbonate fuel cell is a promising device for 
stationary power and heat supply and has been continuously 
developing by many countries including USA, Japan, and Korea. 
It was demonstrated by several manufacturers that in the sub- 
MW region MCFC power plants can reach electrical efficiencies 
of 47%. By making use of the heat generated by the system, total 
efficiencies of more than 80% can be achieved [1]. 

Performance and availability of MCFC stack are greatly 
dependent on its operating temperature. Due to its high operating 
temperature of about 873 K, this concept has certain advantages 
over low-temperature fuel cell systems. The temperature is high 
enough to allow for internal reforming, which means to pro- 
duce hydrogen from different types of fuel gases within the cell 
stack itself. Also due to the high temperature, no expensive cata- 
lysts are required; nickel and nickel oxide are sufficiently active 
as electrode catalysts. However, when the operating tempera- 
ture is higher than 973 K, material corrosion accelerates greatly, 
and electrolyte loss increases, which increases the risk of short- 
circuit and shortens the stack lifespan [2]. Consequently, it is 
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crucial for improving the performance and life of the fuel cell 
that the operating temperature be controlled in an appropriate 
range. 

For getting a better thermal management of MCFC stack, 
it is necessary to develop the thermal model, and find the 
suitable control method based on this model. Over the last 
decade, many researchers all over the world have made great 
progress on MCFC modeling to improve its performance and 
cost competitiveness with energy conversion devices currently 
available. Various mathematical models have been established 
in the research on the internal mechanisms, ranging from a one 
dimensional non-isothermal model to a three dimensional non- 
isothermal and non-isobaric models [3—6], but the most existing 
models of MCFC are not ready to be applied in synthesis. In 
Ref. [7], a nonlinear mathematical model of an internal reform- 
ing MCFC stack is developed for control system. In Ref. [8] a 
control-oriented quantitative relation model of MCFC is devel- 
oped between the temperatures of the stack and flow rates of the 
input gases. However, temperature control is only mentioned 
and the model-based control strategy is not discussed. 

As above mentioned, it is necessary to develop a simple 
model, and find suitable control method for dynamic thermal 
transfer. In this paper, we analyze the dynamic thermal trans- 
fer model of MCFC stack based on the modeling of Shen [8]. 
This model is nonlinear with disturbances and uncertainty as 
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Nomenclature 

A area of heat exchange (m?) 

c mass heat capacity (J kg7! K7!) 

C heat capacity (J K7!) 

F Faraday constant (26.8 A h mol`!) 

j local current density (A m~?) 

k local heat-transfer rate coefficients 
(JST! m? KT!) 

q energy rate available from the input fuel and oxi- 
dant gas (J S7! K7!) 

qmo heat generation of 1 M H20 (J mol`!) 

Q heat energy 

r heat radiation coefficient 

S cross-section area of the flow channel of fuel and 
oxidant gas (m?) 

T temperature (K) 

u flow-rate of gas (m s7!) 

V operating voltage of the cell (V) 

We width of cells (m) 

Greek letter 

p density (kg m7?) 

Subscripts 

e electrode—electrolyte plate 

ef between fuel gas and electrode—electrolyte plate 

eo between oxidant gas and electrode—electrolyte 
plate 

es between electrode—electrolyte plate and separator 

f fuel gas 

inf original quantity of the fuel gas at the inlet 

o oxidant gas 

in_o original quantity of the oxidant gas at the inlet 

S separator 

sf between fuel gas and separator 

so between oxidant gas and separator 


shown later in this paper. However, the performance of gen- 
eralized controller largely depends on the availability of an 
accurate model. As a consequence, the application of con- 
ventional control strategies leads to unacceptable closed loop 
performances. The controller must be robust to uncertainty in a 
non-conservative manner and must meet closed-loop objectives 
such as tracking, regulation and disturbance attenuation. 
Variable structure control systems (VSC) [12,13] constitute 
an important class of control systems. The main advantage of 
VSC is the possibility of obtaining a trajectory describing a new 
type of motion called ‘sliding regime’ or ‘sliding mode’, which 
is not inherent in any of the original structures. Sliding mode 
is defined as the motion of the system state along the switching 
surface that is not a trajectory of any of the structures. In the 
sliding mode, the system dynamics are independent of system 
parameters and external disturbances. As a consequence, based 
on VSC systems some works have been reported for the model- 


following control systems [14-16]. In this paper, the variable 
structure control (VSC) algorithm with a rapid-smooth reaching 
law and rapid-convergent sliding mode [17] is used to control 
the operating temperature within a specified range and reducing 
temperature fluctuation. The input/output linearization method 
is introduced into the affine nonlinear system and then a decen- 
tralized variable structure control system is constructed. Genetic 
algorithm (GA) is used to optimize the parameters of the decen- 
tralized VSC, including the parameters of the sliding mode and 
ones of the reaching law. Satisfactory simulation results are 
obtained. 

The organization of this paper is as follows. The affine non- 
linear model with multi-input and multi-output is given in the 
next section. Section 3 gives the design of the variable structure 
control (VSC) in detail. In Section 4, numerical simulations are 
presented to test the validity of the proposed controller. Finally, 
we present our concluding remarks. 


2. Dynamic thermal transfer model of MCFC stack 
2.1. Description and analysis of MCFC 


A single cell consists of anode, cathode, electrolyte plate, 
separator plates, corrugated plates and current collectors. Fuel 
and oxidant gas flow in the anode and cathode channels along 
the corrugated plates. O2 and CO% in oxidant gas react with 
electron at the cathode, and produce CO32-. CO3? moves 
within the electrolyte plate perpendicularly from cathode side 
to anode side driven by the force of concentration difference. At 
the anode, H3 in fuel gas reacts with CO3°~, which produces 
CO2, H20 and electrons by electro-chemical reaction. Released 
electrons are collected by the current collector and then pass 
the corrugated plate and separator perpendicularly. The top and 
the bottom separators are connected to external load equipment 


[9]: 

Hə + CO3? —> CO; + H20 + 2e7, 
4e~ + 2CO + O2 > 2C03%, 

2H? + O2 > 2H20 


anode : 
cathode : 


overall reaction : 


The temperature of stack is mainly affected by the flow rate 
of inlet gas [11]. Slow flowing of gases leads to an adequate 
reaction, less heat lost and heat produced. On the other hand, fast 
flowing of gases results in an inadequate reaction, much more 
heat flown away with the remnant gases, and much reaction heat 
produced. Therefore, the steady temperature of the stack varies 
with the flow-rate of feeding gas of anode and cathode in a 
complex manner. In general, in order to cool the stack and make 
the fuel to be utilized to the maximum limited fuel utilization 
at most, it is assumed that the oxidant (i.e. compressed air to 
cathode) is abundant. 


2.2. Dynamic thermal transfer of MCFC stack 


In this section, the proposed dynamic thermal transfer model 
of MCFC is based on the following assumptions: 
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Cathode and current collector 


Separator Py Corrugated plate 


Anode and current collector 


Fuel gas 


+— Heat conduction 


——— Heat convective 


----»> Heat radiation 
==> Mass transfer 


Fig. 1. Process of mass transport and heat generation and transfer in a MCFC 
cell. 


e The composition change of fuel and oxidant gas because of 

the shift reactions in the cells can be ignored. 

The stack boundary is isothermal. 

e External reforming, co-flow type stack, and the gases flow in 
smooth-walled channels of constant cross-section. 

e The gases flow fast in short channels, and so the heat conduc- 
tion inside the gases along the flow direction, horizontal heat 
conduction of electrode—electrolyte plates and separators can 
be ignored. 

There are mainly the following heat transfer processes in 
an MCFC stack as shown in Fig. 1: 


(a) the input heat energy of fuel gas; 

(b) heat convection between fuel gas and electrode—electrolyte 
plate; 

(c) heat convection between fuel gas and separator plate; 

(d) the input heat energy of oxidant gas; 

(e) heat convection between oxidant gas and electrode— 
electrolyte plate; 

(f) heat convection between oxidant gas and separator plate; 

(g) heat generation of water in the reaction; 

(h) ohmic heat because of electric current; 

(i) heat radiation between electrode—electrolyte plate and sep- 
arator. 

Mass and energy balance analysis [10] is performed on 

the fuel gas, oxidant gas, electrolyte plate and separator and 
the following relations can be obtained [8]. 


2.2.1. Fuel gas 
The temperature of fuel gas is affected by the heat transfer 
processes (a), (b) and (c). A heat balance is written: 


oT = = Qin f + Qef + Ost (1) 


where 
Oin + = qt(Tin£ — Te), 
Ost = kst A(T; — Tr), 


Qef = kepA(Te — Tp), 


Cf = cpppspw, qf = Cin fPin fSfUf 


2.2.2. Oxidant gas 
The temperature of oxidant gas is affected by the heat transfer 
processes (d), (e) and (f). A heat balance is written: 


G2 = Qin_o + Qeo + Vso (2) 
where 

Qin .o = qo(Tino — To), 
Oso = kso A(Ts — To), 


Qeo = keo A(Te = To), 


Co = CoPoSoW, Jo = Cin_oPin_oSoUl0 


2.2.3. Electrode—electrolyte plate 

The temperature of electrode—electrolyte plate is affected by 
the heat transfer processes (b), (e), (g), (h) and (i). A heat balance 
is written: 


dT, 
Ce a = Ofe + Qoe + Qmo + Qohmic + Ose (3) 


where 


Ofte = ket A(Te — Te), Qoe = keo A(To — Te), 


Qno = =>9m0; Qohmic = JVA, 


OF 


To (as 
s 
Ose = 5.67Fres ( a) = (5) | ‘ Ce = CePeScA 


2.2.4. Separator 

The temperature of electrode—electrolyte plate is affected by 
the heat transfer processes (c), (f), and (i). A heat balance is 
written: 


dT; 
Cs- = Qfs F Qos F Qes (4) 
where 


Ors = ks A(T — Ts), Qos = ks A(T, — Ts), 


Qes = 5.67 Te)" _(%)" C= A 
es = 3.0 /Tes 100 100 , s = CsPsSs 


A number of additional quantities are defined: 


KepA KA keo A kso A 
= ef = sf , ja eo l j= so 
Cr Cr Co Co 
kef A keo A 5.67resA 
A45 = ; 46 = , TAn 
Ce Ce 100°C. 
VA Ago kst A 
Àg = — =—4 Ao = È À = š 
oC, FO = Oo C 
kso A 5.67resA Cin_fPin_£5£ 
ài = , (C= a] = ——_—_, 
Cs 100°C, Cf 
Pp Cin_oPin_oSo 
2 E - 


In MCFC stack, the state variables are fuel gas temperature 
x1, oxidant gas temperature x2, electrolyte plate temperature x3, 
and separator temperature x4, uf and uo are control inputs: 
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e The state variable is 
_ 1 ae T 

X = [x1, x2, x3, x4] = [Tf, To, Te, Ts] 
e The manipulated variable is 

U = [u;, uz)" = (up, uo]" 
e The controlled variable is 

T 
y= [yny] = [x1, xa]" 


Using the above definitions an abbreviated form of equations 
describing the temperature of fuel cell follows: 


: Te x1 
X = f(x) + GWU, y= h(x) = H = | | (5) 


X4 
where 
—(å1 + åà2)x1 + 1x3 + 2x4 
—(A3 + à4)x2 + 3x3 + 4x4 
fœ) = 4 4 
Asx + 6x2 — (Às + à6)x3 + à7xX4 + ABAG 
Aiox1 + à11x2 + à12x3 — (A10 + à11)x4 — A12x4 
a1 (Tin f — x1) 0 
0 02(Tin.o — x 
G@) _ 2( in_o 2) 
0 0 
0 0 


This temperature model will be nonlinear with disturbances 
and uncertainty as shown later. The main reasons are as follows: 
(1) there exists load disturbance (the variation of load current); 
(2) it is difficult to analyze the thermal transfer in a stack due to 
its complicated construction, so accurate thermal transfer model 
may not be established; (3) there are many physical parameters 
to be determined by experiments, and every parameter may be 
different greatly under different conditions. In this paper we 
mainly consider the load disturbance d(x, Aj). 

Variable structure control (VSC) offers insensitivity to 
parameter variations, rejection of external load disturbance and 
fast dynamic response. The temperature model has a normative 
affine nonlinear form that makes it possible to design a variable 
structure controller of the MCFC stack. 


3. Design controller for MCFC 


Our goal is to design decentralized sliding mode control 
laws for nonlinear temperature model described by Eq. (5). We 
introduce the input/output linearization method to the variable 
structure control system, and then a whole decentralized design 
variable structure control system with rapid-smooth reaching 
law and sliding mode is constructed. 


3.1. Linearization 


The input/output linearization method makes the output and 
input in linear relation by state transition. For a complete descrip- 
tion of input/output linearization method and notation used in 
the following expression, see Ref. [18]. 

Consider a kind of MIMO nonlinear system as follows: 


X = f(x) + GWU, y= h(x) 


wherex € R”; u € R”; ye R”; te R; G = [g1 (x) Em(x)] 
is a matrix of n x m; f(x) € R”; g(x) € R” (k=1, 2,...,m). The 
item h(x) is a smooth vector field. 

Now we introduce the symbol of Lie derivative: 


ah aL. h; ƏL! h; 
0 rip — Ot p f kp — Í 
Ppr hpm a p a gt 
i=1,2,...,m (6) 


Then define the structure constant (relative rank) r;: 


mong j-l aL hi 
ri =minş j|Lg Ly hi = Jx g(x) # 0; 
j=1,2,...,n;k=1,2,...,m}, i=1,2,...,m 
(7) 


Get the time differential of yı from Eq. (5): 


ahı (x) 
ox 


j= (f(x) + Gu) = Lip + Lgyuhı (8) 


where 


Lih =[1 0 0 O] fœ 
= fi(x) = —(à1 + à2)x1 + 1x3 + À2x4 (9) 


Levi =[1 0 0 0] Gyu = o1 (Tint — x1)u1 (10) 


Since L4 hı = Lg, Lhi #O(Tin¢ Æ x1), for yı =h; the 
structure constant (relative rank) r1 = 1. 
Get time differential of y2: 


h 
j= RO o + G(x)u) = L hz + LG,ho = Lphħ2 (11) 


where 


Lyh =[0 0 0 1] fŒ = fal) = Atom +à11x2 


+A12x4 — (Ato + à11)x4 — Anox4 (12) 


LG,h2=[0 0 0 1] G@u=0 (13) 
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Sle 


v4 desired final value 


reaching model 


slidind suface 


slidind model 


Fig. 2. Graphical interpretation of VSC. 


Continue to get time differential of y2: 


P aL jhe 1 1 
In = S+G) = LfLpha + Laub phn (14) 
where 
jo i 
aL ¢ho had 
ax 4h7x} , 
—(A1o + à11) — 41213 
2 . 
Lh = LfL}h = á fix) 15) 
i a 47x3 i 


—(à10 + Ani) — 4A 12x3 
Leu pho = Ly Lth + Lg L h 
= Aga (Tint — x1)u1 + A402(Tino — x2)u2 (16) 


Since LGuL ho A 0(Tin £ A X1 and Tino # x2), for y2 =h2 
the structure constant (relative rank) r2 = 2. We define a Jacobian 
matrix J from formulas (10), (13), and (16): 


0 0 
LgiL phi Lg, L phi 
1 1 
Leg, L shz Lg, L pha 


= (17) 
For the actual system (Eq. (5)), J is nonsingular and the charac- 
teristic of the zero-dynamics is globally exponential stable. The 
system can be linearized and stable. 


3.2. Obtain the control law of VSC 


The design problem consists of selecting the parameters of 
each structure and defining the traveling logic. The first step is 
to define a mode s(t) =0, along which the process can slide to its 
desired final value. The structure of the controller is intention- 
ally altered as its state crosses the mode in accordance with a 
prescribed control law. Fig. 2 shows graphical interpretation of 
VSC. 


There are many options to choose the sliding mode; the s(t) 
selected in this work, is a rapid-convergent sliding mode pre- 
sented in Ref. [17]. 

The sliding mode is 


s = ọ(p)e + ce? (18) 


where: g(p)=p+w and p=d/dt is differential operator; y > 0; 
0<fB=81/B2<1, fı and Bo (6; <2) are positive odd num- 
bers and tuning parameter, which helps to define s, and e is 
the tracking error. 

The objective of control is to ensure that the controlled vari- 
able be equal to its reference value at all times, meaning that e 
and its derivatives must be zero. The reaching law [17] is chosen 
as 


$ = —e|s|“sen(s) — ks, e,k > 0,1 >a>0 (19) 


where the —e|s|“sgn(s) item is the sufficient condition for 
s to have sliding mode, so s— 0 in limited time. |s|* acts 
as the “smooth” function and —ks guarantees rapid reaching 
velocity. 

Formulas (18) and (19) have the same structure, with the 
corresponding relation: ko y, æ < B,e< 5. 

Once the sliding mode has been selected, attention must be 
turned to the design of the control law that drives the controlled 
variable to its reference value and satisfies. 

For the MIMO system (5), we define the number of switch 
function equal to the number of control input, m = 2. The switch 
function can be represented as 


s(x) =[s1X) s2] (20) 
The switch surface is 


Si = {x|s(x) =O}, i=1,2 


Using decentralized variable structure control scheme, we 
design the switching function s; according to the structure con- 
stant r; of yj. 

Considering the structure constant r1 = 1, r2 =2 calculated in 
Section 3.1, we get 


sı = pı(p)ei + sief 0 <B<1) (21) 


where 


t 
gi(p) = p! + via, e = 1 (Yı — ya )dt, 
0 


s2 = gr(per + szef (22) 


where t 
gop) = P? + Yap! + yaz, e = f (y2 — ya2)dt 

The yqa=[va1, ya2] is the desired output; Wio=1; Wit, 
Wi2,--+Ji,r; are constants, whose numerical values determine 
the probability of output error on the sliding mode.From 
(5), (17), (21) and (22), we can get the time differential of 
vector s; is 
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5 = A(x) + Ju — p(p)ya 


t p= 
ipl f Q= yapar Qi - ea | 


F 2 
t b-1 
528 | (y2 — vio (y2 — ya2) 


(23) 


where 
[5 mæl 
Nya ply na) 
aw yio=l, y20=1 
Soo, Lha) 
j=0 
Lihi) + yi h(x) 
L4ho(x) + Wo,1L pha(x) + Wo,2ha(x) | 

gi(P)yai 

g(p)ya = 
p2(P)ya2 


Taking (19) into (23), we can get the control law as 


H 2J! i P ee een 


u2 —é|s2|“sgn(s2) — ks2 


t b-1 
spl f (x1 — vara (x1 — ya1) 


t B-1 
528 ff (x4 — vio (x4 — ya2) 


zi —e|s1|%sgn(s1) — ksı 
=J x 
—elsı|*sgn(s2) — ks2 
| LyM + Yih) | 


L7ho(x) + Wo,1L pho(x) + W2,2h2(x) 


pi(P) yal 
p2(p)ya2 


t B-1 
asl f (x1 — yanar (x1 — yaı) 
(24) 


t B-1 
528 ff (x4 — vio (x4 — ya2) 


where L Lh 1(x), Lih(x) and Lh (x) are calculated respec- 
tively by (8), (12), and (15); sı and sz are obtained respectively 
from (21) and (22). 

Here is the input design of the variable structure control of 
nonlinear temperature system (5). The robustness of the nor- 
mative affine nonlinear system with chattering is proved in 
[17]. Selecting parameters of the sliding mode and those of the 
reaching law are considered as one of the main underlying prob- 
lems associated with VSC. This is accomplished by formulating 
the VSC parameters selection as an optimization problem and 
genetic algorithm (GA) are used in the optimization process. 


3.3. Optimize the parameters of controller 


Genetic algorithm (GA) is directed random search techniques 
which can find the global optimal solution in complex multidi- 
mensional search spaces. GA was first proposed by Holland 
[19] and has been applied successfully to many engineering 
and optimization problems [20,21]. It is a searching algorithm 
based on natural selection and species heredity. While simulat- 
ing the hereditary behavior of natural biology, such as aberrance 
and competition, GA improves the solution of a problem, and 
obtains a satisfying or optimum solution. There are basically 
three genetic operators used to generate and explore the neigh- 
borhood of a population and select a new generation. These 
operators are selection, crossover and mutation. More details 
about GA can be found in [22]. 

In this paper, GA is used to optimize the parameters of vari- 
able structure control, including the parameters of the sliding 
mode and those of the reaching law. The parameters searched 
for optimum values are 92,1, 91,2, 61, Ê; 1,1, $2, £, k and a. 
To start the proposed GA method, a performance index must 
be defined. The selection of the performance index depends on 
the objective of the control problem. The procedure described 
the use of GA in optimizing the parameters of variable structure 
control is summarized as follows: 


(1) Determining the ranges of parameters of VSC based on 
experience, and generating the initial population. 

(2) Evaluate some performance index for all optimized param- 
eters generated in step 1. In the present work, the following 
index is selected: 


T 
min f {(x — xaT Wi(x — xa) + u” Wu} dt 
0 


where W, and W2 are the weighting coefficients, and xq is 
the vector of steady state operation. 
(3) Use genetic operators (selection, crossover, mutation) to 
produce new generation of parameters of VSC. 
Evaluate the performance index in step 2 for the new gen- 
eration of parameters of VSC. Stop if there is no more 
improvement in the value of the performance index or if 
certain predetermined number of generation has been used, 


otherwise go to step (3). 
X = f(x) + G(x)U 
X = h(x) 


VO 
linearizing 


(4 


wa 


Fig. 3. Illustration of temperature control. 
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Figs. 4-7. Trajectories of MCFC temperature variable. 


4. Simulation 


In this section, we present numerical experiments to show the 
validation of the proposed variable structure control based on 
the dynamic thermal model of MCFC. The detailed parameters 
about MCFC stack used for simulation can be found in [8]. The 
structure of VSC for MCFC temperature control system is shown 


in Fig. 3. 

In normal working condition, the current density 
jo=1500Am~?. The steady output value of the 
operating temperature of the MCFC stack is xqg= 


[ 900 896 930 921 ] Tk. We assume that there exist 
an initial fuel gas temperature error 20 K Gy = 920K) and a 
step disturbance Aj with the value of 150 A m~? (equal to 10% 
of regular condition) at the fifth minutes. The variable structure 
controller is used to adjust the stack temperature to its steady 
value. 

To start the proposed algorithm, a decision has to be 
made about the GA parameters which include population size, 
crossover probability, mutation probability, and number of gen- 
erations. General guidelines available in the literature can be 


used in the selection process of the GA parameters [22,23]. After 
so many trials, a population size of 90, a crossover probabil- 
ity of 0.7, and a mutation probability of 0.001 are used. Using 
the GA design procedure described in Section 3.3, initial range 
and optimal value of control law parameter that minimize the 
performance indices are, given by 


¢1,1 = 5.6(1-20), sı = 1(0.1-3), 92,1 = 101-15), 
2,2 = 26(0.1-30),  B=0.7(0.1-0.9), c2 = 10.1, 3), 
e = 10.1), k = 6(1-20), «œ = 0.8(0.1-0.9) 


With the control of VSC, we get a group of the MCFC tem- 
perature trajectories. The trajectories of fuel gas temperature 
(x1), oxidant gas temperature (x2), electrolyte plate tempera- 
ture (x3) and separator temperature (x4) are showed respectively 
in Figs. 4-7. Above results show that the temperature can 
track a given temperature reference value well. The temperature 
fluctuations of stack only occur when current density changes 
suddenly, the fluctuations attenuate gradually, and the temper- 
atures achieve the controlling accuracy after a short period of 
time. The upper temperatures are within safety area. The tem- 
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ul (ms!) 


t(min) 


u2 (m s!) 


t(min) 


Figs. 8 and 9. Curves of manipulated variable u; and u2. 


perature error of the fuel gas is less than 11 K, the electrolyte 
plate is less than 1 K and the error of oxidant gas is less than 
12 K, and these conditions are permitted in MCFC stacks. The 
fuel gas temperature and separator temperature x4 (output y2) 
are robust to the current disturbance and the sliding modes sı 
and s2 are invariant to disturbance. From Figs. 8 and 9, veloci- 
ties of fuel gas (u;) and oxidant gas (u2) are within the area of 
control (the fuel gas velocity: 0-2 ms~!, and oxidant gas veloc- 
ity: 0-8 ms~!). Above simulation results show that the control 
system can reach stable finally, and its dynamic characteristics 
are fairly good, and it is robust to the disturbance. 


5. Conclusion 


The operating temperature of the stack is the most important 
variable which is controlled in the MCFC system. However, 
the existing mechanism model is too complicated to meet the 
design requirements of the control system, and less model- 
based control strategy is discussed in past works. According 
to the demand of controlling the temperature of MCFC stack, 
normative affine nonlinear model of MCFC temperature charac- 
teristic considering the disturbance and uncertainty is analyzed 
based on the modeling of Shen [8]. A model-based variable 
structure control (VSC) with GA is presented for MCFC. 
Variable structure control is a special discontinuous control 
technique applicable to a broad variety of practical systems 
due to its easy implementation and good robustness to unmod- 
eled dynamics, parameter variations and external disturbances. 
We introduce the input/output linearization method to the vari- 
able structure control system, and design decentralized variable 
structure control system with optimized parameters by GA. 
Simulation results show that the method presented in the 
paper can control the operating temperature to reach the target 
value quickly, which verified the effectiveness and robustness 
of it. 
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